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It is known that a mesoscopic normal-metal ring pierced by magnetic flux supports persistent current. It was 
proposed recently that the persistent current exists also in a ring made of a band insulator, where it is carried by 
electrons in the fully occupied valence band. In this work, persistent currents in rings made of band insulators 
are analyzed theoretically. We first formulate a recipe which determines the Bloch states of a one-dimensional 
(ID) ring from the Bloch states of an infinite ID crystal created by the periodic repetition of the ring. Using 
the recipe, we derive an expression for the persistent current in a ID ring made of an insulator with an arbitrary 
valence band e{k) . To find an exact result for a specific insulator, we consider a ID ring represented by a periodic 
lattice of TV identical sites with a single value of the on-site energy. If the Bloch states in the ring are expanded 
over a complete set of A*' on-site Wannier functions, the discrete on-site energy splits into the energy band. 
At full filling, the band emulates the valence band of the band insulator and the ring is insulating. It carries a 
persistent current equal to the product of A'^ and the derivative of the on-site energy with respect to the magnetic 
flux. This current is not zero if one takes into account that the on-site Wannier function and consequently the 
on-site energy of each ring site depend on magnetic flux. To derive the current analytically, we expand all TV 
Wannier functions of the ring over the infinite basis of Wannier functions of the constituting infinite ID crystal 
and eventually determine the crystal Wannier functions by a method of localized atomic orbitals. In terms of the 
crystal Wannier functions, the current at full filling arises because the electron in the ring with TV lattice sites is 
allowed to make a single hop from site i to its periodic replicas i ± TV. In terms of the ring Wannier functions, 
the same current is due to the flux dependence of the Wannier functions basis, and the longest allowed hop is 
from i to j ± Int(TV/2) only. Finally, we estimate the persistent current at full filling in a few realistic systems. 
In particular, the rings made of real band insulators (GaAs, Ge, InAs) with a fully filled valence band and empty 
conduction band are studied. The current decays with the ring length exponentially due to the exponential decay 
of the Wannier functions. In spite of that, it can be of measurable size. 

PACS numbers: 73.23.-b, 73.61.Ey 



I. INTRODUCTION 

It is known that a conducting ring pierced by a constant 
magnetic flux can support persistent electron current circulat- 
ing around the ring yj] . The persistent current is well known 
to exist in a superconducting ring, and also in a normal-metal 
ring as long as it is mesoscopic (of size comparable or smaller 
than the electron coherence length). Early work dealing with 
persistent current and flux quantization in superconducting 
rings [2] also contains results relevant to the normal metal 
rings. Later it was mentioned ^ that one can have circulating 
cuiTents for free electrons in ballistic normal-metal rings, and 
the reference lUd proposed that persistent currents should exist 
in the disordered normal-metal rings. The persistent currents 
were indeed observed in the disordered normal-metal rings 
15-9] as well as in ballistic conducting rings 1 10, 11]. Details 
are reviewed in various papers il2l-ll4ll . 

One of us recently proposed fly\ that persistent current ex- 
ists also in a ring made of a band insulator, where it is carried 
exclusively by Bloch electrons in the fully occupied valence 
band. The effect is interesting for two reasons [15J. First, it 
seems to be the only effect manifesting electron transport in a 
fully occupied valence band. In a standard conductance mea- 
surement, the band insulator is biased by two metallic elec- 
trodes and one observes the tunneling of conduction electrons 
from metal to metal through the energy gap of the insulator. 
One does not observe the transport of valence band electrons 



in the insulator Second, in the fully occupied valence band, 
Bloch electrons cannot be scattered by phonons or other elec- 
trons. Therefore, the electron coherence length should be 
huge even at room temperature if the valence band is separated 
from the conduction band by a large energy gap. The persis- 
tent cuiTent in the insulating ring should thus be observable 
up to high temperatures as a temperature-independent effect 
as long as excitations into the (empty) conduction band are 
negligible. By contrast, persistent currents in metallic rings 
are observable only at low temperatures, js^. 

We report on a study of persistent currents in rings made of 
band insulators using the formalism of the Bloch and Wannier 
functions. It was partly motivated by a pioneering analysis 
fl6] of persistent cuiTents in a ID ring-shaped periodic lattice 
composed of the lattice sites with a single value of the on-site 
energy. The authors derived persistent cuiTent / for spinless 
electrons at zero temperature assuming that they move along 
the lattice by means of the nearest-neighbor-site hopping. The 
result reads 
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for even A'^e, where is the electron number in the ring, N is 
the number of the lattice sites, cf) is the magnetic flux piercing 
the ring, (pQ = h/e is the flux quantum, and Fi is the hopping 
amplitude. Equations ([T]i and (|2]l show a nonzero current for 
a conductor {N^ < N). However, they show a zero current 
for an insulator {N^ = N), which contradicts our claim that 
persistent current exists also in an insulating ring. 

Therefore, one of our goals is to give exact proof that a ring 
made of a band insulator does support a non-zero persistent 
current carried by the fully filled valence band. We also pro- 
vide a simple tool allowing to estimate persistent currents in 
rings made of realistic band insulators. We finally present a 
derivation of an exact persistent current expression for a spe- 
cific model insulator. This paper is structured as follows. 

In Sect. II, we formulate a recipe which determines the 
Bloch states of a one-dimensional (ID) ring from the Bloch 
states of an infinite ID crystal created by the periodic repeti- 
tion of the ring. In Sect. Ill we use the recipe to derive an 
expression for the persistent current in a ID ring made of a 
band insulator with an arbitrary valence band e(fc). The ex- 
pression generally shows that the current is not zero, although 
it does turn to zero for a cosine-shaped energy band, which 
was used in the analysis with nearest-neighbor-site hopping 

To derive an exact result for a specific insulator, in Sect. IV 
we consider a ID ring represented by a periodic lattice of N 
identical sites with a single value of the on-site energy. If the 
Bloch states in the ring are expanded over a complete set of 
N on-site Wannier functions, the discrete on-site energy splits 
into the energy band. At full filling, the band emulates the 
valence band of the band insulator and the ring is insulating. 
We find that it carries a persistent current equal to the product 
of N and the derivative of the on-site energy with respect to 
the magnetic flux. An alternative explanation of why Eqs. ([T]i 
and give zero current at full filling is that they do not take 
into account the flux-dependence of the on-site energy. 

To get closer to an exact formula, in Sect. IV we expand 
all N Wannier functions of the ring over the basis of Wannier 
functions of the constituting infinite ID crystal. In terms of 
the crystal Wannier functions, the current at full filling arises 
because the electron in the ring with N lattice sites is allowed 
to make a single hop from site i to its periodic replicas i± N. 
Alternatively, in terms of the ring Wannier functions the same 
current is due to the flux dependence of the Wannier functions 
basis, and the longest allowed hop is from « to i ± Int(7V/2) 
only. 

In Sect. V we eventually derive the crystal Wannier func- 
tions by a method of localized atomic orbitals (LCAO) and 
express the persistent current at full filling by means of an ex- 
act formula. The formula is found to be in good agreement 
with a numerical LCAO calculation. The numerical data are 
provided for a ID ring made of an artificial band insulator, 
namely for a GaAs ring subjected to a periodic potential em- 
ulating the insulating lattice. 

In Sect. VI we estimate persistent currents in the rings made 
of real band insulators (GaAs, Ge, InAs) with a fully filled 



valence band and empty conduction band. The current at full 
filling decays with the ring length exponentially due to the 
exponential decay of the Wannier function tails. In spite of 
that, it can be of measurable size. 

Our results are summarized in Sect. VII. We relate them 
briefly to the paper by Kohn iITtIi . which defines the insulat- 
ing state as a property of the many-body ground state, and 
which contains a remark implicitly admitting the existence of 
nonzero persistent current in insulating rings. 

We note that some of the results of Sect. V were obtained 
in the conference contribution [18] by means of a formally 
different approach. In another conference paper ifioll the per- 
sistent current at full filling was analyzed for a ID ring with 
Kronig-Penney potential. 

II. ELEMENTARY THEORY AND THE GENERAL RECIPE 

To review basic concepts, we consider the circular ID ring 
with circumference L, pierced by magnetic flux 0. The elec- 
trons in such ring are described by the Schrodinger equation 
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ip^ix) = e{(f))ilj^{x), (3) 



where x is the electron position along the ring circumference, 
m is the electron mass, V{x) is an arbitrary potential applied 
along the ring, £(0) is the electron eigen-energy, and ■00 (2;) is 
the electron wave function. Due to the ring geometry, ip^{x) 
has to fulfill the periodic condition 



iIj4,{x) = i)^{x + L). 
We define the wave function ^P(j,{x) by transformation 
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ip^{x) = exp i 
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Setting ip^{x) = eyiY){—i^-^x)ip^{x) into ([3) and dUi we 



obtain the Schrodinger equation 



with the boundary condition 



ip^{x) = e{(t>)(p^{x) 



(6) 



ip^{x + L) = exp{i2TT—)(p^{x). 
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Here the magnetic flux enters only the boundary condition. 
In the ring geometry, V{x) obeys the periodic condition 



Vix) ^ V{x + L), 
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which is the same as for the infinite ID crystal with lattice 
constant L. Therefore, equation ^ with condition dHJ is 
mathematically identical with the Schrodinger equation 
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where V{x) is the infinite periodic potential with period L, 
obtained by the periodic repetition of the ring potential V{x). 
Equation (|9]l has the well known Bloch solution 

(pk{x) = exp{ikx)uk{x), (10) 

where the function Uk (x) fulfills the periodic condition 

Uk{x) ^ Uk{x + L), (11) 

and k is the electron wave vector from the interval (—00,00). 
Clearly, the wave function Q and Bloch solution ( fTOl i coin- 
cide for k = ^4-- Let us discuss this coincidence in detail. 

Consider first the ring with zero magnetic flux. To obtain 
the wave functions in such ring, it suffices to take the Bloch 
function ( fTOb and to restrict it by the periodic condition 



(pkix) = (pk{x + L), 



(12) 



which is the condition (|7]l for (j) — Q. Due to the condition 
(fTSI i. the wave vector k becomes discrete: 
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Thus, in the ring with zero magnetic flux and specified poten- 
tial V{x), the eigen-function tfn{x) and eigen-energy e„ can 
be calculated simply by setting k — into the Bloch so- 
lutions (pk{x) a e{k), calculated for the same potential V{x) 
repeated with period L from x = —00 to x = 00. This recipe 
can be generalized to nonzero magnetic flux as follows. 
Arbitrary magnetic flux cf) can be written in the form 



where (p' is the reduced flux from the range < 
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alternatively from < 0, (j)o), and n is one of the values n — 
0, ±1, ±2, . . .. Setting (Hill into Q one can write (0) in the 
form 
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where the function ipnff,.{^) 
periodic condition ijj'^ ^, (x) 
condition d?) now reads 



exp {i^nx') -00 (2;) obeys the 
<h- + L). The boundary 



Vn.d>. (x + L) = cxp(i27r— , (x). (16) 
fpo 

Similarly, in the Bloch function theory it is customary to 
express the wave vector k by means of the relation 



k = — ri + k , 



(17) 



where k is the reduced wave vector from the first Brillouin 
zone and the integer n plays the role of the energy band num- 
ber. Using ( [TT] ) we can write the Bloch function ( fTOb as 



Vn,k. {x) = exp (ifc'a;) u'^_^, (x), 



(18) 



where the function ^, (x) = exp (i^nx'j Uk{x) fulfills the 

periodic condition u ,, (x) — u , (x + L). It is then easy to 
verify that 

fnM. (x + L) ^ cxp{it L)ip'^ f^, (x). (19) 

The equations ( fTSl l and (fT6] l coincide with the equations 
(fTsT l and ( fT9] l, respectively, if 

L 00 ' 
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One can thus formulate the following general recipe. In 
the ring with a known potential V{x), the eigen-function 

(p'^ ,j,.{x) = (p^{x) and eigen-energy e„(0') = can be 
calculated by setting ( l20b or ( |2TI ) into the Bloch solutions 
(fin.k-ix) = fk{x) ae„(fc') = £{k), calculated for the same 
potential V{x) repeated with period L from x — —00 to 

X = 00. 

The general recipe holds for any potential V{x) obeying the 
cyclic condition It therefore holds also for any potential 
which additionally obeys the periodic condition 

V{x)^V{x + a), (22) 

where a = L/N and N is the (integer) number of periods 
a within the period L. This type of potential will be con- 
sidered in the rest of the paper. Figure [T] illustrates how an 
infinite ID crystal is created from a ID ring with lattice pe- 
riod a and length L = Na. Such crystal is still described by 
Schrodinger equation (|9]l, except that now V{x) — V{x + a) 
and consequently Ufc (x) = Uk{x + a). However, the condition 
Uk{x) = Ufc(a; + L) holds as well because F(x) — V{x + L). 
We can rederive the general recipe briefly, if we compare 
Uk{x) = Uk{x + L) with ip^{x) = ip4,{x + L) and ipk{x) = 

exp{ikx)uk{x) with ip^{x) = exp (i^-^xj ip4,{x). We see 

that (p^{x) coincides with <pfc(x) for k = 

It remains to express the persistent current. The Bloch elec- 
tron in state (n, k ) moves with velocity w„(A: ) = 7i~~a^^- 



We set into the Bloch electron velocity the formula (I20b . We 
obtain the expression w„(0 ) = '^^a^^ ^ ' which is the elec- 
tron velocity in state {n, </> ) in the ring. The current carried 
by the electron in state (n, (j) ) reads 



/n(0 ) = - 



evn{(j>') den{4>) 



L d(j)' ' 

The total persistent current circulating in the ring is 



d 



(23) 



(24) 



where we sum (at zero temperature) over all occupied states 
n = 0, ±1, ±2, ... up to the Fermi level. In the following text 
we use the formulas (l23T l and ( |24] | with symbol 4> changed to 
(f>, where (f> e< -^,^) or alternatively (j) e< 0,(f>o)- 
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FIG. 1. The top figure shows schematically the ID ring composed of 
the periodic lattice with A*' lattice sites (A'' = 3), periodicity a, and 
length L — Na. For simplicity, the lattice potential of the ring is 
modeled by a square-well potential V{x). The bottom figure depicts 
the infinite ID periodic potential obtained by the periodic repetition 
of the ring potential. 



III. PERSISTENT CURRENT IN ID RING MADE OF A 
CRYSTAL WITH ARBITRARY VALENCE BAND e{k) 

Consider first an infinite ID crystal created from a ID ring 
with the lattice period a and length L — Na as it is illustrated 
in Fig. [T] Let ey{k) be the energy dispersion of the valence 
band of Bloch electrons in that crystal. The wave vector k in 
e„(fc) is supposed to be from the Brillouin zone < >■ 
In what follows we skip the valence band index (v) for sim- 
plicity and we use notation e(fc). We can express e{k) by 
means of the infinite Fourier expansion 

e(fc) — ao + ai cos(fca) + 

a2 cos(2fca) + . . . ajv cos{Nka) + . . . , 



where 



and 
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Now we apply the general recipe from Sect. II. We set for 
k in Eq. the equation fc = + ^). We obtain the 

expression for the ring energies en(0): 



£„((/)) = e{kn{(f>)) = ^ ftj COS (fc„(0) ja) 

3=0 



where 
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As shown in appendix A, expansion (|28] l can also be written 
in the form 



M 



£n{4') — COS {kn{<p) ja) + ilj sin (fc„(0) 

(30) 



where M = Int(iV/2) and the Fourier coefficients are given 
by equations ( 11021 ) and ( 11031 ). Expressions ( |28] | and ( l30b are 
mathematically equivalent, but expression ( |30l ) reveals an im- 
portant difference from the Bloch energy ( |25] ). Generally, the 
Bloch energy dZSl l consists of an infinite number of Fourier 
terms, with the coefficients aj being independent on mag- 
netic flux (p. On the other hand, expression ( l30b shows that 
the Bloch energy in the ring can be expressed through a fi- 
nite number of Fourier terms, with the Fourier coefficients il 
depending on (f> [see Eqs. ( 11021 ) and ( I103H . 

In the rest of this section we use Eq. ( 1281 ). The same results 
can be obtained by means of Eq. dSOl l. but derivations are 
more cumbersome. 

Setting Eq. ( |28] ) into Eq. ( l24b we obtain the persistent cur- 
rent 
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If the ring contains A^e spinless electrons and Ne is odd, 
summation over n in Eq. dSTI) includes the occupied states 
n = 0, ±1, ±2, . . . , ±(A''e — 1)/2. After simple manipulations 
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Performing the summation over n we arrive to the formula 
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j=N,2N,... ^ 
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where G< ^). Note that the first sum in ([33]l omits 

the terms j = N, 2N, 3N, . . . , included in the second sum. 

If Np is even and G< 0, (po), we sum in Eq. ( [3TT l over the 
occupied states n = 0, ±1, ±2, . . . , ±(iVe/2 - 1), -Ne/2. 
Similar manipulations as before lead to the result 
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where </> €< 0, ^o). The results ^ and ^ hold for the 
conducting (A^e < A*") as well as insulating {N^ — N) rings. 

For = N the first sum on the right hand side of Eqs. 
(33[ and ( |34l i becomes zero and both equations take the form 



j=N,2N,... ^ 



(35) 



where ip G< —^,^) and N can be even as well as odd. 
Equation ( |35] | expresses the persistent current in a ID ring 
made of a band insulator with an arbitrary valence band Sk- 

An important property of the formula ( l35T l can be recog- 
nized without specifying the dispersion law eu- The disper- 
sion e{k) is a periodic function of period 27r/a. As long as it is 
analytic in the whole Brillouin zone < — - - >, the Fourier 
coefficients aj in the expansion ( l25T l obey the inequality ll20ll 

Ifljl < Ce-^\ (36) 

where C and ^ are positive constants. Equation (|35] | combined 
with inequality (|36] | shows that the persistent current in the 
insulating ring is in general not zero albeit it decreases with 
the ring length at least exponentially. In the next sections, 
we will demonstrate for various specific insulating rings the 
exponential decrease, and we will see that the current can be 
of measurable size in realistic systems. 

The pioneer results ([T]) and (|2]) can be obtained from Eqs. 
( |33] | and (|34] | if we keep only the term j — \ and we put 
fli = — 2ri. Equations ([T]i and Q predict for insulating rings 
(iVe = N) a zero persistent current, which contradicts equa- 
tion (|35] |. Origin of this contradiction is evident. Keeping 
only the term j ~ I means that only the nearest-neighbor- 
site hopping is operative. This approach, also known as tight- 
binding approximation, simplifies the general dispersion law 
(l25l l to the cosine shaped dispersion e{k) oc cos(fca). How- 
ever, equation (l35T l shows that nonzero persistent current in 
the insulating ring of length N is (mainly) due to the Fourier 
coefficient which is in general not zero for a realistic dis- 
persion law. We will see in the next section that physically 
represents the hopping amplitude for a single hop from site i 
to its periodic replicas izL N. 



IV. PERSISTENT CURRENT IN INSULATING ID RING IN 
THE WANNIER FUNCTIONS FORMALISM 

We consider again the ID ring with lattice period a and 
length L = Na, sketched in Fig. [T] Simultaneously, we con- 
sider the infinite ID crystal obtained by the periodic repetition 
of the ID ring, as it is illustrated in Fig. [T] For the sake of 
clarity we recall a few major points. First, the Bloch func- 
tion tfk{x) and Bloch energy eik) in the infinite crystal are 
described by Schrodinger equation 



Hipk{x) = e{k)ipk{x) , 



H 



2m dx 



Second, according to the general recipe of Sect. II, the Bloch 
function and Bloch energy in the ring, Lpn,iji{x) and £„(</)), are 
given by equations 

H^nA^) = <P/c=fe„(0)(a;), en(0) = e{k = kn{(f>)) , 



(38) 



Third, the Bloch vector k is assumed to be from the Brillioun 
zone < — -f,^ >, or alternatively from < 0,^ >. This 
means that e{k) is the dispersion of a single energy band and 
ifk {x) are the Bloch functions of that band. As before, the 
functions ipk{x) and e{k) do not involve any band index, and 
e{k) is supposed to represent the valence band of the crystal. 
Below we define the Wannier functions by restricting us to the 
Bloch states of that valence band. 



A. Wannier functions of the infinite crystal 

We first define the Wannier functions for the crystal. Ac- 
cording to a standard textbook definition, the Wannier func- 
tion W{x — la) at the lattice site x = la is defined as 



W{x - la) = Wiix) = ^ I e-^'='°$fc(a;)dfc , (39) 

ZTT 



where ^k{x) is the Bloch function of the infinite crystal. We 
note that ^k{x) is the same Bloch function as ipk{x) in equa- 
tion (IJTI i. The only difference is that ^k{x) is normalized as 



/ I 'S>k{x) P dx^l 
Jo 



(40) 



while the normalization of ipk (x) in the infinite crystal is un- 
specified. In the following text, represents the nor- 
malized Bloch function of the infinite crystal and the symbol 
tfk{x) with k = kn{4>) represents the Bloch function of the 
ring. [In the ring (pk{x) is normalized by equation (l48Tl.1 

Since the crystal is infinite, fc is a continuous variable and 
the Bloch functions fulfill the orthogonality condition 



dx = —Sik' - k) . (41) 



One can then easy verify the orthogonality condition for the 
Wannier functions. 



^+V{x), V{x)^V{x + a). (37) 



Wi*,{x)Wi{x)dx^Sij, . (42) 

Finally, definition ( |39] | can also be written in the inverted form 

^k{x)=J2 e'''"Wi{x) (43) 

where one sums over all lattice points of the (infinite) crystal. 



6 



We define the matrix elements Tj by equations 



Ti'.z = - / Wi.{x)HWi{x)dx ^ Ti^i, , r, = r,_,,, 

J — oo 

(44) 

where /, /' — 0, ±1, ±2, • • • ± oo and it is taken into account 
that Ti'j depends only on the difference j = I — I'. The Bloch 
energy e(fc) can be written as e{k) = {^k\H\^k) / {^k\^k) ■ 
If we set for $fe equation (|43] |. we easy find 



Fj cosjka. 



(45) 



Comparison of Eq. i45[ with Eq. (l25l l shows that aj = 2rj. 
Obviously, 

= I ' <k) cos{jka)dk, (46) 

7r/2a 

which is an alternative form of equation (l44l l. 

Finally, setting aj = into the result (35i we get 



j=N,2N,... ^ 



(47) 



Equation ( l47T i together with Eq. (|44] | express the persistent 
current in the insulating ring through the Wannier functions of 
the constituting infinite crystal. We defer discussion of these 
expressions to subsection C, where equation (|47] ) is derived 
from the Wannier functions of the ring. 



B. Wannier functions of tlie ring 

Let us define the Wannier functions of the ring. In the 
ring with N lattice sites and a fixed value of magnetic flux 
only the discrete Bloch vectors fc„(0) are allowed, where 
kn{0) e< -f,f > or alternatively fc„(0) e< 0,^ >. In 
both cases there are N allowed values of n (specified below), 
which means that the originally continuous band e(fc) consists 
of N discrete levels £„((/)). The eigen-functions of these lev- 
els represent a complete set of N Bloch functions obeying the 
orthogonality condition 



Na 



(pl,(x)(fkix) dx = 6k,k' = Sn 



(48) 



where k = kn {4>)- Consequently, in the ID ring with A'^ lattice 
sites the Wannier function at site x = la can be defined as 

wix - la) = wi{x) = ^Yl ' (49) 



where k — kn{(l>) and summation over k means the summa- 
tion over the N allowed values of n. Definition (|49] | coincides 
with definition (3.7) in Kohn's paper IitIi . except that Kohn 

assumes k — A:„(0). 



In the following calculations we choose the Brillouin zone 
fc„(0) G< 0, ^ >, for which the allowed n are simply 



n^O,l,...,N -1 



(50) 



without regards to the parity of N and sign of 

Using Eqs. (|48] | and (|49] | we can easy verify the orthogo- 
nality condition 



Na 



wv{x)*wi{x) dx = 5v^i 



(51) 



Finally, definition ( |49] l can also be written in the inverted form 



1 

(^) = ^ E • (52) 



In this paper the Wannier functions wi {x) are often called the 
ring Wannier functions while the Wannier functions Wi [x) are 
called the crystal Wannier functions. The ring Wannier func- 
tions are magnetic-flux dependent, the crystal Wannier func- 
tions are not. 

We note that the results obtained in the rest of this section 
can also be obtained by choosing /c„ (0) e< — f , f >■ In that 
case the allowed values of n are 

n = 0,±l,±2,...,±(A^-l)/2 (53) 

for odd N, while for even N one has 

n = 0,±l,±2,...,±(7V/2-l),-iV/2 for > (54) 

and 

n = 0,±l,±2,...,±(iV/2-l),+A^/2 for (/)< 0. (55) 

Evidently, it is more straightforward to work with Eq. ( |50l l. 

Now we express the persistent current at full filling in terms 
of the ring Wannier functions. We write £„(0) in the form 



Na 



en[(l)) = / Lpl ,j,-.{x)Hipk^(^^){x)d: 



(56) 



We need to rewrite the right hand side of Eq. ( |56] | in terms of 
the matrix elements 



T, 



I'l 



Na 



Wi, (x) H wi {x) dx , 



(57) 



where 1,1' ~ . . . N —1. We set into Eq. ( |57] | the equation 
gUl. We find 



v. I 



iV EE 



^ikn' {<p)l' a ^ — ikji{(p)la 



No 



^ / 'PL'{^){^)Hfk„w{x) dx 



N-l 



ra=0 



(58) 
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Since Ti'j depends only on the difference between I and I', 
we can write 



Ti,,i = Ti,_i = T, = le*^^^' e^^"-'£„(0) . (59) 
Inverting Eq. ( l59b we get 



N-l 



N-1 
i=0 



(60) 



By means of Eq. ( l60l l. the ground state energy of the fully 
filled valence band is readily obtained in the form 



Af-l 



= ^£„(0)-iv-ro(</)) 



(61) 



Finally, the persistent current at full-filling reads 



dE 
"d4 



-N 



dTo 



(62) 



One can see from Eq. dSTl ) that the matrix element To rep- 
resents the on-site energy of the Wannier state wi {x). Since 
wi{x) depends on magnetic flux, so does Tq. This however 
means that the current ( l62b is in general not zero. We have 
thus arrived at the same conclusion as in the preceding sec- 
tion, where the persistent current in the insulating ring (Eq. 
[35] l was derived from the Bloch states and general recipe. 



C. Ring Wannier functions in tlie basis of crystal Wannier 
functions 



The result ( |47] |. or equivalently the result ( l35T l. can be ob- 
tained from the result ( |62] i. if the ring Wannier functions are 
expanded over the basis of the crystal Wannier functions. To 
show this we proceed as follows. 

The Bloch functions of the ring are given by Eq. ( l38b and 
normalized as (p1{x)ipk{x) dx = 1, where k = kni4>) 
[see Eq. (|48]l1. The Bloch functions 4>k{x) coincide with 
ipk {x) except for the normalization constant. As they are nor- 
malized by means of Eq. ( |40l l. in the ring we have the relation 



(t>k{x) , 



(63) 



We set for in Eq. (|63] l the definition ( |43T l. Then we set 

Eq. ( |63] l into the right hand side of Eq. ( |49] l. We obtain 



1 



^'(^) = ^E E e-*'=('-'')'^W^p(a:) , (64) 

k l' — — oo 

where k = and summation over k means the summa- 

tion over n = 0, 1, . . . , A^— 1. Using substitution I' — Nr + s, 



where s = 0, 1, . . . (A^ - 1) and r = 0, ±1, ±2, • • • ± oo, we 
further obtain 



Af-l 



r— — OO s— 
N-l 

-T 

N ^ 



iV 



n=0 



Eventually 



(65) 



(66) 



Equation (l66l l expresses the ring Wannier function at the ring 
site I, wi{x), through the crystal Wannier functions Wi'{x). 
For = equation ( |66] | coincides with equation (5.9) in 
Kohn's paper 1 17]. 

Setting Eq. (|66] | into Eq. ( fSTl ) we obtain 

Tj = / w* {x) H Wo {x) dx 



oo oo 



-(r-r') 



Na 



WNr'+jix)HWNr{x)dx 



r— — oo r — — oo 

oo , oo rNa-rNa 



^Y e'^^^^o^ 

S— — 00 r= — oo ' 



0-rNa 



WNs+ji^)HWo{x)dx 



s—~-oo 



Y e I WNs+j{x)HWo{x)dx 



(67) 

By means of the matrix elements Tj, equation (|67] ) can be 
written as 



(68) 



s— — oo 

For j = equation ( [68] t gives 

oo , 

To = To + 2 y TsN cos(2^^s) . (69) 

We substitute Eq. f6% into the right hand side of Eq. ( |62l ). 
We obtain the persistent current at full filling in the form 



Attn 



Y, * ^sN sin(27r— s) 



(70) 



or equivalently in the form ( |47] | which coincides with the re- 
sult ([Mil for = 2rj. 

In the rest of this section we discuss the physical meaning of 
the coefficients Tn, r2N, etc., appearing in expression ( l47l) . 
Since aj ~ 2rj, the meaning of the coefficients a^, a2N, etc., 
is the same and does not need an extra discussion. 
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According to the definition (|44] |. in the infinite ID crystal 
the matrix element F^v represents the hopping amplitude of 
a single hop from the Wannier state at site I to the Wannier 
state at site I + N. In the ring with N lattice sites the site 
I + N is a periodic replica of site /, which means the same 
site. Appearance of Fat in expression ( |47] | therefore implies, 
that the electron in the ring with N sites can make a single 
hop from the Wannier state at site / to its periodic replica at 
site I + N. In other words, the electron hops from site I to the 
same site by making one round around the ring, and Fat is the 
hopping amplitude of the process. Note that we speak about 
hopping in the ring, but in terms of the Wannier states of the 
infinite crystal (c.f. Eq. l44l i. 

Similarly, the coefficient F2Ar describes the process in 
which the electron hops from site I to the same site by mak- 
ing two rounds around the ring. It follows from relation ( |36] l 
and equation aj = 2Fj, that F2Ar is much smaller than Fjv- 
Its contribution to the right hand side of Eq. ( |47] | is therefore 
negligible. The meaning of the coefficients F3JV, F4JV, etc., is 
analogous and their contribution is negligible even more. 

We can summarize the result ( |47T i as follows. The persistent 
current in the insulating ring with N lattice sites, derived in 
terms of the crystal Wannier functions, is due to the fact that 
the electron at the ring site I is allowed to make a single hop 
from site I to its periodic replicas I ± N. (The sign minus is 
allowed as well due to the symmetry reason.) 

The process in which an electron at the ring site I hops from 
site / to the same site by making one round around the ring 
may seem meaningless, because such process does not exist 
in the basis of the ring Wannier functions. In that basis the 
energy spectrum £„(</)) is given by equation ( I6OI 1. As shown 
in appendix B, equation (|60] l can be rewritten for odd N as 



' 



(71) 



and for even N as 



E 



-ikn{<t>)ja 



(72) 



where T^j = T* . According to the definition (|57] |. the matrix 
element Tj represents the amplitude of hopping from the ring 
Wannier state wi{x) to the ring Wannier state wi'{x), with 
j ^ I — V being the hopping length. The maximum hopping 
length between the ring Wannier states, allowed by equations 
(ItTI i and (|72] |. is j = ±Int(Af/2) only, and the current in 
the insulating ring (Eq. |62] | is due to the fact that the ring 
Wannier states depend on magnetic flux. Hopping from I to 
I ± N, I zL 2N, etc., emerges after the ring Wannier states are 
expressed (see Eq. l64l i through the flux-independent crystal 
Wannier states. This provides a different picture of the same 
physics. 




2d 



FIG. 2. Left: One-dimensional model of a single isolated atom, de- 
scribed by equation J73b . The atomic potential, v{x), is modeled 
by the potential well of width 2d, embedded between two infinitely- 
thick potential barriers of height Vo- Right: The potential V{x) in a 
ring-shaped periodic lattice composed of N one-dimensional atoms 
(A*' — 6). The lattice period is a, the ring circumference is L = Na. 
In the numerical calculations of Sect. V, the parameters of the ring 
are chosen to emulate the GaAs ring with conduction electrons sub- 
jected to the periodic potential of A'^ quantum dots: the electron ef- 
fective mass is m = 0.067mo, the period of the quantum-dot lattice 
is a = 40nm, the ring length is L = Na, and the size and depth of 
the potential well (quantum dot) are 2d = 15nm and Vo = 18meV, 
respectively. For these parameters the potential well supports one 
bound state, the ground state with energy Sa ~ — 10.4meV. The 
ground-state orbital ipaix) is sketched schematically. 



V. THE LCAO APPROXIMATION, APPLICATION TO THE 
ID LATTICE WITH RECTANGULAR POTENTIAL WELLS 



To evaluate the persistent current (Eq. |47j for a specific 
crystal lattice, we need to determine the crystal Wannier func- 
tions of the lattice. As an example we consider the ring- 
shaped ID lattice in Fig. |2] The lattice is composed of the 
artificial ID atoms with a single atomic level ea- The atomic 
orbital (pa{x) and energy ea in the isolated atom positioned at 
X — obey the Schrodinger equation 



2m dx'^ 



v{x) 



(Pa{x) = eaiPa{x) , 



(73) 



where v{x) is the atomic potential modeled as a rectangular 
potential well centered at a; = (figure |2] left). 

In accord with the general recipe of Sect. II, assume that 
the ring with N atoms (figure|2] right) is periodically repeated 
as shown in figure [1] This repetition gives rise to the infinite 
ID crystal with potential V{x) which is periodic with period 
a and ranges from x = —00 to x = 00. Evidently, V{x) is a 
sum of all isolated atomic potentials, that means 



V{x) = 



/ — — oc 



v{x — la) 



(74) 



where la is the position of the l-th lattice point. The Bloch 
function ^k{x) and Bloch energy e(fc) of the electron mov- 
ing in such periodic potential can be expressed analytically in 
the approximation of localized atomic orbitals (LCAO). The 
Bloch function is approximated by the LCAO ansatz 112 III 



(75) 



I— — 00 
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where Lpix — la) is the atomic orbital at the lattice point la 
and C is the normalization constant (given below). The Bloch 
energy e{k) can be expressed as 



(76) 



where 



ft2 d 



^ = -^^ + ^(a:), Vix) ^ v{x~ la). (77) 



1— — 00 



If we use the LCAO ansatz (fTSb . a simple textbook calculation 
1 211] gives the formula 



where 



7o + 2 V"!, 7, cos (kja) 
1 + 2 2^j^i aj cos [kja) 



dxipa [x - ja)ipa (x) , (79) 



7j = - / dxipaix-ja)V'{x)ipaix), 7o=7j=o, (80) 
and 



= V{x) ~ v{x), 



(81) 



with w (a;) being the atomic potential at the lattice site ^ = 0. In 
Appendix C we express aj and jj analytically for the model 
atomic lattice of Fig. |2] 

Using the general recipe of section II, we obtain the spec- 
trum of the ring, en{<l>), directly from the Bloch energy dTSl) . 
It reads 



en(0) 



70 + 2 T,JLi 7i cos [kn i(t))ja] 
1 + 2 J2T=i [knWja] ' 



(82) 



where A;„ (0) = + ^ ) ■ One can substitute the spectrum 

(|82| i into the persistent current formula (|24] | and one can eval- 
uate the sum on the right hand side of Eq. ( l24l i numerically. 
We will refer to this approach as to the numerical LCAO ap- 
proach. The approach provides a straightforward numerical 
evaluation of persistent current in insulating rings [N^ — N) 
as well as in conducting rings (N^ < N). It does not provide 
understanding of the effect, however, it is useful for verifica- 
tion of analytical results. 

We mainly want to verify our result for insulating rings, the 
formula (|47] i. For a meaningful comparison with the numeri- 
cal LCAO approach, the hopping amplitude Tj in the formula 
(|47| i has to be determined in the LCAO approximation. 

We start by specifying the constant C in the ansatz ( fTST l. 
Setting Eq. (fTsT l into the normalization condition (|40] | we find 



C 



2 aj cos jka = N{k) 



(83) 



The same result is obtained if we set Eq. ( |75] ) into the orthog- 
onality condition (|4TI ). 

The crystal Wannier functions are defined by equation ( [39l ). 
We set into equation the ansatz (|75l) with C = 1 I jNjk ). 
We obtain the crystal Wannier function in the form 1221 1: 



Wi{x) ^ X! ci,_i ipa{x - I'a) 



where 



a 

2^ 



cos (kna) 



dk 



(84) 



(85) 



Expression (|84] | is the crystal Wannier function in the LCAO 
approximation. It is easy to verify that the expression (|84| ) 
fulfills the orthogonality condition W\\ exactly. 

If we express the crystal Wannier functions in Eq. ( l44b by 
means of equation ( |84] |. we find that 



E 



E 



(86) 



Equation (I86t is the hopping amplitude Tj in the LCAO ap- 
proximation. An alternative form of Eq. ( |86] ) is the equation 



a 

2^ 



r/a 

— n/a 

70 ■ 



dk cos{jka) 

2 E/^i 7j' cos (fc/g) 
1 + 2 X] "i' cos (fc/a) 



(87) 



obtained by setting the Bloch energy jTSl l into the Fourrier 
transformation ( |46] l. 

We also note that the Bloch energy jTSl ) is an alternative 
form of the expression ( |45T l with coefficients Fj given by 
equation ( |86] l, since both forms follow from the same LCAO 
approximation. Proof of coincidence of both forms is straight- 
forward, but lengthy and therefore not shown. 

The LCAO expression (l86b together with expressions (|85] | 
and dSSl l are still complicated for analytical calculations. For- 
tunately, the expressions can be simplified, if we take into ac- 
count that the coefficients aj and jj decay with increasing j 
exponentially. Exponential decay is due to the fact that the 
atomic orbitals (pa are exponentially localized at the lattice 
sites, which is the case for any realistic orbital. If the local- 
ization is strong enough, then 1 ^| ai |^| a2 |^ • ■ • and 

1 »| 7i |>| 72 |> In that limit the expressions ( |85l ) 

and ( |86] l are calculated in appendix D. The results are 



c„ ~ (-ai)" \An\, n > , c_ 
where Aq = 1 and 



■n — : 



(88) 
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and 

~ -7i(-ai)(j-i) , j>l. (90) 

Expressing Tj in equation ( |47] l by means of Eq. ( |90] l and 
keeping only the term with j = N we obtain the result 

/(^) = - ^ 71 ^ sin f2^^') . (91) 

The last equation is the persistent current in the insulating ring 
in the LCAO approximation. The following features of the 
result ( |9T| l are worth stressing. 

For positive ai equation ( |9TT i predicts the exponentially de- 
caying current with alternating sign, I{N) ^ {—ai)^^^^\ 
Just this result holds for the ring in figure |2l Indeed, ai is 
positive for a positive (pa (x) such like the ground state orbital 
in figure 121 which follows from definition ( |79] l as well as from 
the final expression for aj (appendix C). 

For negative ai the dependence I{N) ^ {—aiY^^^^ im- 
plies the exponentially decaying current with fixed sign. For 
the square-well model in figure|2]the coefficient ai is negative 
for the first excited bound state (not shown in the figure), be- 
cause the atomic orbital of the state is an antisymmetric func- 
tion of X with respect to the well center. Thus, if the atomic 
level Ea and atomic orbital (pa{x) belong to the first excited 
bound state, the sign of the current is the same for any N. 

Equations ([T]) and (|2]i show that the persistent current in the 
conducting ring {N^ < N) exhibits the alternating sign in 
dependence on the parity of N^. This result is independent 
on the specific properties of the conduction band. If the ring 
is insulating (N^ — N), the result (|9T| i shows that the sign 
of the current with respect to N is either alternating or fixed, 
depending on the properties of the valence band. More pre- 
cisely, this statement holds for the model ring in figure |2l We 
will see in the next section that the sign behavior can be even 
more complicated in rings made of the real band insulators. 

Comparison of equation ( |9T] l with Eqs. ^ and (|2]l also 
shows another important difference between the insulating 
and conducting rings. In the former case the current decays 
with N exponentially, in the latter case it decays like 

Origin of all these differences is easy to see. Equations 
([T]i and (|2]i depend only on the hopping amplitude Fi, while 
the result is determined by the hopping amplitude Fjv. 
The amplitude F n involves overlap of the Wannier functions 
Wi{x) and Wi+n{x). If we express the coefficient q/_; in 
equation ( |84] | by means of equations dSST l and ( |89] |, we obtain 

oo 

Wi{x) ~ I ^K'-'l I (-"i)"'"" I M^-l'a) ■ (92) 

^ — — 00 

If we keep x — na and I = for simplicity, we obtain 

Wo{na)c^^l={-ai)\--\^a[0), \n\>l. (93) 
^/tt I n I 

It can be seen that the Wannier function decays with distance 
from site / = exponentially [24] and even shows the alter- 
nating sign for positive ai. All these properties are reflected 
by persistent current at full filling. 




FIG. 3. Persistent current in the insulating ID ring as a function of 
magnetic flux for various ring lengths N . The ring is the periodic lat- 
tice with single-level potential wells (Fig. |2]l which contains Ne 
spinless electrons. It is insulating due to full filling {Ne = A*'). The 
full lines are the results of the LCAO formula \9\\ and the dashed 
lines are the results of the numerical LCAO approach. The ring pa- 
rameters used in these calculations are specified in the caption of 
figure [2l The coefficients otj and 7j in Eqs. ( 19 U and ( |82t are given 
by expressions presented in appendix C. 



Finally, we calculate the persistent current quantitatively for 
the numerical parameters in figure [D These parameters em- 
ulate the ID GaAs ring modulated by a periodic lattice of N 
quantum dots. In such lattice, the discrete energy level of the 
quantum dot splits into the energy band. Therefore, the ring 
behaves at full filling (iVg = N) like a band insulator and at 
partial filling (N^ < N) as a metal. 

The figure [3] shows the persistent current in the insulating 
ring as a function of magnetic flux for various N. The results 
evaluated by means of the LCAO formula (1911 1 are compared 
with the results obtained by the numerical LCAO approach. 
The numerical LCAO approach nicely confirms basic features 
of the formula ( |9T1 ) including the alternating sign of the cur- 
rent. The quantitative difference between the numerical and 
analytical LCAO data could be suppressed by choosing the 
ring with a more strongly localized atomic orbitals. 

The dependence on the ring size is presented in figure |4] 
The left panel shows the persistent current as a function of 
in the insulating ring and in the conducting ring at half filling. 
The conducting ring shows for large N the decay / cx 1/N, 
in the insulating ring the cuiTent decays with N exponentially. 
For example, the insulating ring with N = 5 supports the per- 
sistent current as small as ^ O.OOlnA. Such small persistent 
cuiTents are detectable by sensitive techniques |9, 25.. 26] . 

The right panel of figure |4] shows the transition from the 
conducting state to the insulating state at = N, achieved 
by varying for fixed N. Experimentally, could be var- 
ied say by using a suitably designed metallic gate, while the 
quantum-dot lattice might be realized either by means of the 
periodic array of gates or by producing the ID ring with a 
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FIG. 4. Persistent current in the same ring as in the preceding figure 
for magnetic flux (j> = 0.2500 • The left panel shows the dependence 
on A'^ for the insulating ring (Ne — N) and conducting ring at half 
filling (Ne = N/2). The right panel shows the dependence on Ne for 
various values of A*'. In both panels, the data for N,; = N shown by 
open symbols were obtained by means of the LCAO expression J91t 
and the data for partial filling (A^e < N) were obtained by means 
of expressions (T) and The full circles in the left panel show 
the data for N,; = A'^, obtained using the numerical LCAO approach. 
Note that only the size of the current, | / |, is shown: for all presented 
data the sign of the current varies with varying parity of Ne . 



periodically alternating cross-section size. 



VI. ESTIMATES OF PERSISTENT CURRENT IN RINGS 
MADE OF REAL BAND INSULATORS 



In the previous section, we have discussed the ring made of 
the artificial band insulator, namely the ID GaAs ring with a 
few conduction electrons subjected to the periodic quantum- 
dot potential. Now we want to estimate persistent current in 
rings made of the real crystalline band insulators (with fully 
occupied valence band and empty conduction band). The 
problem is that such rings are three-dimensional and the ID 
theory from the preceding text is not applicable directly. 

The problem is outlined in figure|5] The top right sketch in 
the figure is a 2D sketch of the 3D ring (or a hollow 3D cylin- 
der) of width w, created from the cubic 3D lattice of atoms. 
Along such ring there is no periodicity with the lattice period 
a, there is only periodicity with L. However, the theory in 
Sects. Ill - V was developed for the ID rings with periodicity 
a along the ring circumference. The question is how to extend 
the ID theory to be applicable to the ring-shaped 3D cluster 
of atoms sketched in the figure |5] 

In principle, for L ^ it is reasonable to write the 3D 




FIG. 5. The top figure shows the 3D ring (on the right hand side) 
created from the cubic 3D atomic lattice (on the left hand side) with 
the lattice period a. The ring of width w is defined by two concentric 
circles with the center positioned at random. The resulting ring is a 
3D cluster of atoms which can no longer be viewed as the ID lattice 
with period d. The bottom figure shows the ring artificially linearized 
by means of elastic deformation. Repeating the linearized ring with 
period L in the x direction and with period w in the y direction, we 
obtain the infinite crystal. In the limit L w the effect of the elastic 
deformation has to disappear and the Bloch states in the obtained 
infinite crystal have to coincide with the electronic states in the real 
ring (top right sketch). This allows us to use the equations l |94| ) and 
([95}. We call this model the model of the linearized 3D cluster. 



version of equation ( l24l i in the form 

d 



del) 



(94) 



with the spectrum of the ring, e„(0, ky = 0, k^), determined 
as follows: 

£7i{<P, ky = 0,fc^) = e{kx = ^{n+-^),ky = 0, fc^), (95) 

L (po 

where £{kx, ky, kz) is the Bloch energy of the 3D crystal cre- 
ated by periodic repetition of the linearized 3D cluster (Fig|5]i. 
Note that the unit cell of the crystal is the linearized cluster 
The Bloch energy £{kx, ky, k^) should therefore not be con- 
fused with the Bloch energy of the crystal with unit cell a^. 
We call this approach the model of the linearized 3D cluster. 

The model also allows to rewrite for 3D rings all equations 
of Sect. III. First of all, since the period a along the ring no 
longer exists, it has to be replaced by the period L (replace a 
by L and by unity). Second, the Bloch energy e{kx) has 
to be replaced by the Bloch energy £{kx, ky — 0, kz) just in- 
troduced. Third, the current ( l35T l has to be summed over the 
occupied states kz. Clearly, the model of the linearized 3D 
cluster allows to avoid a full 3D treatment of the ring-shaped 
3D cluster. However, the model is still cumbersome for prac- 
tical use because the Bloch spectrum £{kx, ky, kz) needs to be 
calculated for a nontrivial artificial crystal resulting from the 
periodic repetition of the linearized cluster 

A more simple approach is to consider the artificial 3D ring 
obtained by elastic bending of the bulk crystal, as it is shown 
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FIG. 6. The 3D ring of length L and width w, created by elastic 
bending of the cubic 3D atomic lattice with lattice period a. In the 
limit w <^ L the effect of the elastic bending has to disappear and the 
electronic states in such ring have to coincide with the Bloch states 
in the constituting 3D crystal as described in the text (equation |97}. 
We call this model the model of the elastically bent 3D crystal. Obvi- 
ously, the resulting 3D ring shows for w; ^ L an artificial periodicity 
with period a, which is not the case for the realistic 3D ring (fig. |5}. 



in figure |6] In this case we have instead of the formula (|94] l 
the formula 



d 



(96) 



where the energy spectrum ky, kz) is given by equation 



£n{_4'j ky^kz) — £{kx — {n H ■ ) , ky, k^ 
I\ a 00 



(97) 



with e{kx, ky, kz) being the Bloch energy in the bulk crystal 
before it is bent. We call this approach the model of the elas- 
tically bent 3D crystal. In this model there is an artificial peri- 
odicity with period a along the ring. In spite of this, the model 
can provide a reasonable quantitative information on how the 
persistent current in the 3D ring depends on the realistic band 
structure of the constituting bulk crystal. In the model, the per- 
sistent current at full filling is given by equations (|47] | and (|46] | 
modified as follows. The Bloch energy £{kx) in the Fourrier 
transformation ( |46] | has to be replaced by the Bloch energy 
^(fcx, fcy, kz) and the current Wh should be summed over the 
occupied states ky and kz ■ Practical calculations are relatively 
straightforward, because the Bloch spectrum of the bulk crys- 
tal, e{kx^ ky, kz), is known for many insulators. 

In what follows we apply the model of the elastically 
bent 3D crystal to the rings with thickness equal to the unit 
cell of the bulk crystal. The energy dispersion of the ring, 
£n{4i,ky — 0,kz — 0), is extracted from the bulk crystal 
dispersio n e( kx, ky = 0,kz ~ 0), calculated by microscopic 
methods 11281 12911 . The obtained persistent current should be 
viewed as a minimum value estimate, because the thickness 
of realistic rings is much lar ger t han the unit cell. 

Using standard methods i28l Eoll . in the figure [T] we nu- 
merically reproduce the well-known band structure of the 
band insulators like the GaAs, Ge, and InAs. Specifically, 
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FIG. 7. The panels in the right column show the valence bands 
of various band insulators like the GaAs, Ge, and InAs. Specif- 
ically, the e(fc) dispersion curves of the light-hole band (full 
line), double-degenerate heavy-hole band (dashed line), and low- 
est band (dotted line) are plotted for the fc-values on the line 
(L[—TT/a, —ir/a, — 7r/a], r[0, 0, 0], L[7r/a, n/a, ir/a]). The panels 
in the left column show the hopping amplitude Fjv, calculated sep- 
arately for each band: the e{k) dependence of the band is Fourier- 
transformed by using the equation i46l . The full, dashed, a dotted 
lines show Fjv for the light-hole, heavy-hole, and lowest bands, re- 
spectively. When compared with Fjv of the light-hole and heavy- 
hole band , Fjv of the lowest band has opposite sign because of the 
opposite curvature of the e(fc) curve and the presented data should be 
multiplied by —1 as shown in the figure (but see the comment 1.27.1 ). 



the figure summarizes the efc curves of the light-hole band, 
heavy-hole band and lowest band, calculated for the line 
{L[— n/a,— n/a,— n/a], T[0, 0,0], L[n/a, n/a, n/a]) in the 
first Brillouin zone. We set each of these e{k) curves numer- 
ically into the Fourier transformation (l46T l and we calculate 
for each band the hopping amplitude Fjv- The resulting Fjv 
dependencies are shown in the figure |7] 

By means of the figure|7]we estimate the persistent current 
as follows. We form the ring by elastic bending of the crystal. 
We choose the ring orientation for which the k states on the 
line {L[—n/a, —n/a, —n/a],T[0, 0, 0], L[7r/a, n/a, n/a]) 
are directed along the ring circumference. If the ring is 
one-dimensional (with thickness equal to a single unit cell 
of the crystal), the persistent current is given by the ID 
formula (|47] | with the hopping amplitude Fjv given by 
numerical values in the figure |7] As expected, F^r decays 
with increasing N exponentially, however, for the light-hole 
band the decay is much slower than for other two bands. 
Therefore it is sufficient to calculate the persistent current in 
the light-hole band. The results are shown in the figure |7] 
We recall that these results represent the persistent current 
carried by electrons in the fully occupied valence band. The 
following features are worth to point out. 

The persistent current in the insulating ring decreases with 
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VII. SUMMARY AND CONCLUDING REMARKS 
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FIG. 8. Persistent current versus the ring length A'^ in the ring made 
of the real crystalline band insulator (GaAs, Ge, In As). The persis- 
tent current is estimated by means of the ID formula i47l . where Tm 
is the hopping amplitude of the electron in the light-hole band, cal- 
culated in the preceding figure. The lattice constants of the GaAs, 
Ge, and InAs crystals are ~ 0.6nm, i.e., the maximum ring lengths 
considered in the figure are L = Na ~ 120nm. The rings of such 
length could be realizable by means of advanced nanotechnology. 



the increasing ring length exponentially, but a careful choice 
of the insulating material allows to achieve the persistent cur- 
rent of measurable value for a technologically realizable ring 
size. Indeed, the maximum ring lengths considered in the fig- 
ure [8] are close to the value ~ 120nm, in principle achiev- 
able by modern nanotechnology. Concerning the amplitude 
of the persistent current, the value ~ O.OOOSnA estimated for 
the InAs ring of length ^ 120nm is very small, but in prin- 
ciple detectable say by the experimental technique of Refs. 
SEHES- The considered insulators are ordered as GaAs, 
Ge and InAs for the persistent currents ordered increasingly. 
This suggests a simple rule: the smaller the light-hole effec- 
tive mass the larger the persistent cuiTent. 

We conclude with a few remarks showing that the above ID 
estimates hold roughly also for the 3D rings. Similar calcu- 
lations as in the figure [T] were performed for the k values on 
the line (X[-27r/a, 0, 0], r[0, 0, 0], X[27r/a, 0, 0]). We have 
found the energies e(fc) and amplitudes Tn (not shown) very 
similar to those in the figure |7] with the resulting persistent 
currents being of the same order of magnitude as those in the 
figure [8] Moreover, we have observed similar results when 
testing the k values on the lines {K, F, K) and (W, F, W). The 
model of the elastically bent crystal thus provides roughly the 
same persistent cuiTent independently on the fact which crys- 
talline axis is chosen to be directed along the ring circum- 
ference. This suggests that anisotropy of the valence band 
with respect to the k direction does affect the resulting persis- 
tent cuiTent dramatically. Therefore, the estimated values of 
the current should hold reasonably also for the realistic ring- 
shaped crystals similar to the ring sketched in figure |5] 



We have analyzed theoretically persistent currents in rings 
made of band insulators. We started by formulating a recipe 
which determines the Bloch states of a one-dimensional (ID) 
ring from the Bloch states of an infinite ID crystal created by 
the periodic repetition of the ring. Using the recipe, we have 
derived an expression for the persistent current in a ID ring 
made of an insulator with an arbitrary valence band e(fc). 

To find an exact result for a specific insulator, we have con- 
sidered a ID ring represented by a periodic lattice of N identi- 
cal sites with a single value of the on-site energy. If the Bloch 
states in the ring are expanded over a complete set of N on- 
site Wannier functions, the discrete on-site energy splits into 
the energy band. At full filling, the band emulates the valence 
band of the band insulator and the ring is insulating. 

We have found that the ring carries a persistent current 
equal to the product of N and the derivative of the on-site 
energy with respect to the magnetic flux. This current is not 
zero because the on-site Wannier function and consequently 
the on-site energy of each ring site depend on magnetic flux. 

Further, we have expanded all N Wannier functions of the 
ring over the basis of Wannier functions of the constituting 
infinite ID crystal. In terms of the crystal Wannier functions, 
the cuiTent at full filling arises because the electron in the ring 
with N lattice sites is allowed to make a single hop from site 
i to its periodic replicas i zL N. Alternatively, in terms of 
the ring Wannier functions the same current is due to the flux 
dependence of the Wannier functions basis, and the longest 
allowed hop is from ito i± Int(A^/2) only. 

We have eventually derived the crystal Wannier functions in 
the LCAO approximation, and have expressed the persistent 
current at full filling by means of an exact formula. The valid- 
ity of the formula was verified by comparison with the numer- 
ical LCAO approach. We have presented quantitative results 
for a ID ring made of an artificial band insulator, namely for 
a GaAs ring subjected to a periodic potential emulating the 
insulating lattice. 

Moreover, we have provided estimates for the rings made 
of real band insulators (GaAs, Ge, InAs) with a fully filled 
valence band and empty conduction band. The current at full 
filling decays with the ring length exponentially due to the 
exponential decay of the Wannier function tails. In spite of 
that, it can be of measurable size. 

We have ignored the effect of nonzero temperature. We 
expect the persistent currents in rings made of the band in- 
sulators to be temperature-independent as long the thermal 
energy fc^T is much smaller than the energy gap separating 
the valence band from the conduction band. For example, the 
energy gap in the InAs crystal is ^ 400meV. Hence, the zero- 
temperature values of the persistent current, predicted for the 
insulating InAs ring, could persist even at room temperature. 

We have also ignored the electron-electron interaction. The 
many -body studies llSol, 31 1 suggest that in the ID ring consid- 



ered in this paper the repulsive electron-electron interaction 
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would only renormalize the electron transmission through the 
potential barrier separating two potential wells. That means 
that the electron-electron interaction would not essentially 
change our results. 

Finally, we relate our results to the paper by Kohn flTll. 
which defines the insulating state as a property of the many- 
body ground state (see also the review by Resta ll32ll ). Kohn's 
paper contains a discussion implicitly admitting the existence 
of a negligibly small persistent current in the insulating ring 
[see the discussion of equations (2.23) - (2.31) in that work]. 
Kohn points out the exponential localization of the electron 
states in the ring. However, he clearly states by writing his 
equation (2.30), that the current is zero only approximately. 
Just this negligibly small current is calculated in our paper 
Surprisingly, it can be of measurable size. 
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APPENDIX A: DERIVATION OF EQUATION <30) FROM 
EQUATION (28) 



We start by writing Eq. ( I28l l as 

oo 

= ^ a-,' cos {kn{(f>) j'a) , (98) 

j'=0 

Using substitution j' = Ns + j, where s = 0,1,2,... and 
j = 0, 1, . . . (A^ - 1), we obtain 

N-l oo 

= ^ ^asjv+jC0s[fc„(</))(s7V+j)a]. (99) 
Assuming that N is odd we rewrite Eq. ( |99] l as 

oo 

en(^) = X! "'"^ '^'^^ [kn{(p)sNa] 

M oo 

+ ^^a^w+jcos [kn{(l)){sN + j)a] 

j=l s=0 
N-l oo 

+ ^ ^asN+j cos[kn{(p){sN + j)a] , 

j=M+l 5-0 

(100) 



where M — Int(A^/2). Equation dlOOI t can be rewritten as 

OG 

£n{<l>) = ^ asN COS [kn{(l))sNa] 

s=0 
M oo 

+ X! X^l^'^^+J [kni4>)isN + j)a] 

j = l 5=0 

+ a(^sN+N-j) COS [fc„((/)) (sN + N -j)a]j . 
UtiUzing definition ( |29] l we find from Eq. jlOll ) the result 

M 

(101) 

where 

^oi.4>) = y^asN cos(27r— s) , 

oo ^ 

f^f (0) = «j + X! cos(27r— s) [asN+] + clsN-]] , 



s=l 



^^Acf)) = sin(27r— s) [agN-j - asN+j] 



(102) 



We have sofar assumed odd A^. For even N, derivation is 
similar and we find again the results dlOll ) and ( 11021 ). except 
that for j = N/2 we have instead of Eqs. (11021) the equations 



^w/2(<A) = °^sJV+Ar/2 cos(27r— s) , 

s=0 

oo ^ 

^w/2('/') = - as7v+JV/2sin(27r— s) 
Equation (llOll l coincides with Eq. ( |30] l. 



(103) 



APPENDIX B: VARIOUS FORMS OF EQUATION (60) 

We first assume odd N. We can rewrite Eq. ( |60) as 

M N-l 
3 = 1 3=M+1 

(104) 

where M = Int(N/2). Further, 

M M 
M 

(105) 
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where we have used the definition (|29]l and relations 



In addition, for j = 



T* ^ T^j , Tn-j = e"" *^T* , (106) 

which follow from Eq. ( |59] l for real £„(</)). For even N we 
obtain in a similar way the equation 

Af-l 
J = l 

(107) 

Equations d 1051 1 and ( 11071 ) can be written in the unified form 

M 



£„(</>) = Ao + ^|Af cos(fc„(<?!)) ja) + Ajsin(fc„(0) ja)| , 

(108) 



where 



Ao = To , 

Af = (T, + r;), 

A^ = 1(^,-77) 



(109) 



for odd as well as even N, except that for even N and 
j — N/2 we have instead of equations ( |109t the equations 

(110) 



Equation ( IIO8I 1 is formally equivalent with Eq. (llOll l. Coin- 
cidence of both equations becomes obvious if we realize that 
equations ( 1109b and ( lUOI l coincide with Eqs. ( 11021 ) and ( 11031) . 
To observe this coincidence, it is sufficient to express Tj in 
equations ( 11091 ) and ( II 101 ) by means of expression (l68T l. and 
to use the equations ag = To and = 2Tj. 



APPENDIX C: ANALYTICAL FORMULAS FOR AND 7^ 

For the square- well potential in figure [T] the atomic orbitals 
(fia {x — la) can be found easy and the integrals ( |79l ) and dSOl l 
can be calculated analytically. For j = 1, 2, . . . we obtain 



1 + Bd 



S(ja - 2d) 



(111) 



1 + Bd 



-2e, 



« + 2dB{j - 1) 



where 



A = 



1 2m 
IP" 



(M) + ea), 



2rn 
"TP" 



(112) 
(113) 



COs2(Ad) 2B(a-2d) l__e 

7o = Vq— 



-4Bd 



1 + Bd 



-2Ba 



(114) 



Expressions (II 1 lb and (II 12b are not exact (we have skipped 
some exponentially small terms), but they are in excellent 
agreement with numerical integration of the integrals ( |79l ) and 
(ISOl l. As expected, the coefficients aj and jj decay with in- 
creasing j exponentially. The simplest forms are 



7j =i 71 e 



-U~l)Ba 



ai 



ai e 



j>l. (115) 



APPENDIX D: APPROXIMATE FORMULAS FOR 
COEFFICIENTS c„ AND Fj 

Equation dSST l defines c„ by means of the integral which can 
be calculated analytically in a certain limit. In what follows 
the calculation is presented for positive n because c„ = c_„. 
We first express the normalization constant 



/)v(fc) 



1 + 2^ Qfj cosjfca 



E 

s=0 



where 



A, 



(-1)^ 



2 aj cos jka 



2s 



(116) 



(117) 



A closer inspection of the right hand side of Eq. ( II 161 ) shows, 
that the equation (II 161 ) can be approximated as 



v/]v(fc) 



s=0 



As [2ai cos kaY 



(118) 



if aj falls with j exponentially as in equation ( II 15b . and if 
a? ^ 02 ■ We note that the latter condition takes the form 
ai ^ exp(— _Ba) for aj derived in the preceding appendix. 
In general, both conditions are fulfilled simultaneously if the 
localization of the atomic orbitals is strong enough in compar- 
ison with the lattice constant a. 

We set the expansion (II 181 ) into the formula (|85] | and we 
perform the integration. We get 



ail 



s=0 



where 



It can be seen that 



cos(nx) cos'^{x)dx . 



0, 



(119) 



(120) 



(121) 



16 



if s < n or if (s — n) is odd. It can also be seen that 



2 ^n + 2l 



(122) 



if s = n + 21, where ^ = 0, 1, . . . . 

We set equations ( 1122b and ( I121l i into the equation dl 191 ) 
We get an approximate expression for c„. 



n + 2/ 



The last equation can be further rewritten as 



(123) 



i=0 

oo 



n + 21^ 



1=0 



I 



(124) 



where the right hand side is a properly chosen upper limit (be- 
cause An+2l/An < 1). 

If we use the relation 



1=0 



21 



1 



I J y/l -Aal \ l + ^l-4a^ 



the inequality ( I124I I reads 



Cn < (-ai)" \An\ 



(125) 



- 4af V 1 + Vl-4a? 



(126) 



For ai <^ 1/2 it can be further approximated as 

c„<(-ai)"|A„|(l + a?)" (127) 



We return again to the equation ( I123l l and we replace its right 
hand side by a properly chosen lower limit: 



> i-a^r\A,,\{l + aiy' , 



where we have utilized the inequality 

An+2i (n^ 2l\ ( n 



An V / 



> 



(128) 



(129) 



Notice now that the right hand sides of the inequalities ( 1127b 
and ( 1128b coincide. This implies that 

Cn = (-ai)" (1 + a?)" \An\ ^ \An\ . (130) 

Now we attempt to find an approximate expression for the 
coefficient Fj. Equation ( l86b can be rewritten by means of 
substitution n' — n = m as 

oo oo 

= E (eattm - 7m) ^ CnCj+m+n ■ (131) 
m— — oo Ti— — oo 



Since c„ = c_„ , a„ = a_„, and 7„ = 7_„, we obtain 

oo 

= {saao - 7o) cqCj + ^ c„ (cj+„ + Cj_„ 

oo 

+ E (^attm - 7m) [cq (Cj+„ + Cj_„) 
m— 1 

oo 



(132) 

Owing to the exponential dependence (11301) . the last equation 
can be simplified as 



~ E {SaOLra - 7m) ^ C 



■n^j—m — n • 



(133) 



m=0 



Setting for c„ and An the equations (1130b and ( 1117b we have 



Tj ~ E (ea^m - 7m) (-"l) 



j-m 



(134) 



Expressing 7„ and by means of the equation dl 151 ) we fur- 
ther obtain 

{sa - 7o) (-ai)^ 
+ ^(eaai-7i)e-(™-i)^''(-ai)^-™ . 

m— 1 

(135) 

If Q!i ^ exp(— i?a), we can skip all terms with m > 2. Thus 

~ (-ai)^-' [-71 + 7oai)] ~ - (-ai)^"' 71 • (136) 
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